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SUMMARY 


I . 

This  report  describes  the  operation  and  research  activities 
at  the  Norwegian  Seismic  Array  (NORSAR)  for  the  period  from 
/ \ October  1977  to  31  March  1978. 

The  performance  of  the  NORSAR  Detection  Processor  has  been 
somewhat  improved  relative  to  the  previous  reporting  period 
(uptime  increased  from  88.8%  to  91.6%),  although  the  Special 
Processing  System  (SPS)  has  caused  a considerable  number  of 
breaks,  the  longest  one  lasting  for  nearly  5 days.  The  regular 
operation  of  the  new  Event  Processor  (AUTOEP)  was  resumed 
as  of  1 October  1977,  and  the  results  are  now  published 
in  a NORSAR  Monthly  Bulletin  which  is  distributed  to  about 
60  recipients.  Statistics  from  the  first  half  year  of  onera- 
tion  show  an  average  of  12.3  reported  events  per  day,  which 
is  63%  of  the  number  reported  prior  to  the  array  reduction 
The  operation  of  the  data  center  is  working  fairly  well  with 
one  manned  shift.  The  maintenance  contract  with  IBM  was  re- 
duced by  1 October  1977,  and  the  NORSAR  personnel  have  thereby 
increased  their  engagement  in  problems  tied  to  malfunctioning 
of  the  IBM  equipment  (tape  drives,  SPS,  etc.).  The  performance 
of  the  data  communications  systems  (including  the  ARPA  network) 
can  be  characterized  as  good,  although  the  number  of  outages 
increased  during  the  last  3 months  of  the  period.  The  Detection 
Processor  has  not  been  subject  to  major  changes,  however, 
the  AUTOEP  system  has  been  significantly  improved  in  order 
to  meet  the  various  user  requirements.  The  performance  of 
the  array  instrumentation  has  been  stable  and  satisfactory, 
and  the  main  channel  characteristics  show  very  little  change 
from  previous  periods. 

The  research  activities  are  described  in  11  separate  subsections 
of  the  last  chapter  of  this  report  and  cover  research  conducted 
under  NTNF ' s contract  with  ARPA  as  well  as  research  projects 
sponsored  by  Norwegian  authorities.  The  first  subsection  is  on 
the  work  of  the  seismological  expert  group  established  by 
the  Conference  of  the  Committee  on  Disarmament  (the  CCD)  of 
the  United  Nations.  The  second  one  deals  with  statistical 
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models  for  seismic  magnitude,  and  the  following  two  are  on 
P-wave  amplitude  anomalies  and  inversion  of  travel  time 
data.  Then  follow  two  contributions  presenting  results  from 
the  new  microearthquake  array  in  Svalbard,  one  deals  with 
microearthquake  surveillance  and  one  with  the  teleseismic 
detectability  of  the  array.  The  next  subsection  describes  a 
maximum  likelihood  method  for  epicenter  location  based  on  S- 
P time  differences.  Then  follow  one  report  on  precisely 
located  earthquakes  in  the  vicinity  of  the  NORSAR  array,  and 
one  on  macroseismic  data  collection  using  newspaper  ads.  The 
next  one  deals  with  the  seismicity  of  East  Africa,  and  the 
last  subsection  discusses  lithosphere  thickness  in  the 
general  NORSAR  siting  area. 

H.  Gjrtystdal 


OPERATION  OF  ALL  SYSTEMS 


II . 

II. 1 Detection  Processor  (DP)  Operation 

There  have  been  111  breaks  in  the  otherwise  continuous  ooera- 
tion  of  the  DP  system  in  this  reporting  interval.  The  uptime 
percentage  is  91.6%,  as  compared  to  38.8%  for  the  last 
reporting  period  (April-September  1977).  vig.  n.i.i  and  the 
accompanying  Table  II. 1.1  both  show  the  daily  DP  downtime 
for  the  days  between  1 October  1977  and  30  March  1973.  The 
monthly  recording  times  and  up  percentages  are  given  in 
Table  II.  1.2.  As  can  be  seen  from  Table  II. 1.1  the  dominant 
component  governing  the  DP  system  performance  is  the  SPS.  Of 
the  111  breaks  in  this  period,  70  were  caused  by  this  unit 
alone.  Also,  three  of  these  breaks  lasted  for  more  than  two 
continuous  days,  with  the  longest,  117  hours,  break  from 
24  February  to  1 March.  The  breaks  can  be  grouoed  in  the 


following  categories: 

a)  SPS  malfunctioning  70 

b)  Maintenance  stops  9 

c)  Error  on  the  Multiplexor  Channel  3 

d)  Stops  related  to  system  operation  8 

e)  Hardware  problems  5 

f)  Power  jumps  and  breaks  4 

g)  Stops  related  to  program  changes  4 

or  tests 


h)  Stops  related  to  possible  program  2 

errors 

i)  Magnetic  tape  drive  problems  1 

Apart  from  Category  a) , the  numbers  in  the  other  categories 
are  relatively  normal.  The  high  number  of  maintenance  stoos 
(9)  is  partly  related  to  a CPU  error  on  the  360-A  machine, 
but  also  reflects  the  extra  effort  from  NORSAR  oersonnel 
doing  preventive  maintenance  on  the  SPS  unit. 

The  total  downtime  for  this  period  was  363  hours  If)  minutes. 

The  mean-time-between-f ailures  was  1.5  days,  which  is  the 
same  as  for  the  earlier  reporting  period  (April-Sentember  1977) . 


D.  Rieber-Mohn 
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Month 

DP  Uptime 
(Hrs ) 

DP  Uptime 
(%) 

No.  of 

DP  Breaks 

No.  of  Days 
with  Breaks 

DP  MTBF* 
(Days) 

Oct 

728.1 

97.9 

20 

14 

1.4 

Nov 

657.1 

91.3 

16 

14 

1.7 

Dec 

690.2 

92.8 

26 

16 

1.1 

Jan 

640.0 

86.0 

30 

17 

0.9 

Feb 

559.9 

83.3 

10 

9 

2.3 

Mar 

729.4 

98.0 

9 

9 

3.4 

Total 

Period 

4004.8 

91.6 

111 

79 

1.5 

* MEAN -TIME- BETWEEN- FAILURES  = (Total  uptime/No.  of  Up  Intervals) . 


TABLE  II.  1.2 

» 

Online  System  Performance  t ^ 

October.  1977  - March  1978 

j 

J 

1 

II. 2 Event  Processor  Operation 

4 

The  regular  operation  of  the  Event  Processor,  using  the  AUTOEP, 
was  resumed  as  of  1 October  1977.  The  results  are  now  published 
in  a NORSAR  Monthly  Bulletin,  usually  issued  within  two  weeks 
after  the  last  data  date,  and  distributed  to  about  60  recipients. 

Some  statistics  from  the  first  6 months  of  operation  are  given  in 
Table  II. 2.1,  where  it  is  seen  that  12.3  events  have  been 
reported  every  day  in  average.  This  is  62%  of  the  number  of 
eve*.  , reported  during  the  same  6 months  in  1975/76,  a drop 
which  mainly  reflects  the  reduction  of  the  array  from  22  to 
7 subarrays.  When  more  data  are  available,  we  will  look  more 
closely  at  this  drop  in  detectability,  as  well  as  the  location 
accuracy,  and  compare  with  expected  results. 


P.  Engebretsen 
H . Bungum 
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Table  11.2.1 


1 

Tol obo i smte 

Core 

Sum 

Dai  1 y 

) Oct 

77 

215 

165 

380 

12.3 

Nov 

77 

162 

52 

214 

7.) 

Doc 

77 

301 

32 

335 

10.8 

Jan 

78 

156 

54 

210 

6.8 

Fob 

78 

140 

114 

26  3 

9.4 

Mar 

78 

798 

57 

855 

27.6 

II. 3 N0RSA9  Data  Processing  Center  (NDPC)  Operation 
Data  Center 

The  operation  of  the  data  center  is  still  working  fairly  wel 1 
with  just  one  manned  shift,  although  routine  jobs  cover  more 
than  1/4  of  the  shift.  The  users  have,  however,  learned  to 
operate  the  computer  so  that  they  can  run  their  jobs  outside 
the  manned  shift,  if  necessary. 

The  DP  uptime  for  the  period  is  91.6%  and  although  it  is  better 
than  the  last  half  year,  the  number  of  stoos  and  breakdowns 
of  the  SPS  have  increased.  There  have  been  three  major  break- 
downs on  the  SPS  and  those  breakdowns  stand  for  6.8%  of  the 
downtime,  the  other  SPS  stops  1%  and  other  reasons  9.6%.  The 
number  of  stops  outside  office  hours  caused  by  the  SPS  is  ^7. 
This  is  an  increment  of  more  than  50%  compared  to  the  last 
half  year. 
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Since  the  maintenance  contract  with  IBM  was  reduced  (by 
nearly  50%)  1 October  1977,  NORSAP.  personnel  have  solved 
most  problems  in  connection  with  tape  units.  Also  the  EOC 
equipment  has  been  taken  care  of.  In  most  problems  with  the 
IBM  360/40  (B)  and  the  SPS  the  same  personnel  have  been 
involved,  and  to  a certain  degree  reduced  the  IBM  engagement. 

J.  Torstveit 
O . A . Hansen 

Data  Communications  (National) 

The  first  3 months  were  characterized  by  relativelv  few  outages, 
both  with  respect  to  group  as  well  as  single  circuits.  In 
the  last  quarter  the  number  of  outages  increased  for  both 
categories.  Simultaneous  outages  for  groups  of  circuits  will 
almost  always  be  caused  by  carrier  frequency  equipment. 
Approximately  60  outages  were  observed  in  the  reporting 
period,  of  which  March  alone  had  around  40.  Just  a few  outages 
of  this  kind  exceed  1 hour  in  duration.  The  large  majority 
are  observed  over  one  or  two  16-minute  intervals. 

Single  subarray  communications  circuits  have  also  been  affected 
by  the  usual  reasons  such  as:  cable  damages,  intermittent 
equalizer/amplifier  operation?  level  fluctuations,  etc. 

Subarrays  particularly  affected: 


02B 

week 

12, 

5.0% 

02C 

2, 

10.9% 

04C 

If 

2, 

6.6% 

04C 

II 

7, 

24.2% 

04C 

If 

13, 

20.0%. 

I 
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In  addition  the  above-mentioned  and  the  remainina  subarrays 
have  been  subject  to  outages  several  times,  but  the  figures 
are  lower  and  may  vary  between  1 and  3%. 

Modems  and  related  equipment  are  still  in  good  condition.  Just 
a few  times  it  has  been  necessary  to  replace  defective  cards 
(separation  filters).  New  equalizers/amplifiers  were  planned 
by  NTA  to  be  installed  in  November  1977.  Due  to  delays  in 
delivery  the  equipment  is  still  not  installed. 

Table  II. 4.1  shows  outages/degraded  performance  with  respect 
to  subarray  communications  circuits. 


II. 4 The  ARPA  Subnetwork  (TIP  to  TIP,  i.e.,  TIP  incl.  modems, 
lines  and  interfaces) 

The  London  Communications  Circuit 

Apart  from  line  level  adjustment  (10/24)  and  loss  of  carrier 
(11/09),  reliable  performance.  The  Norwegian  Telegraph  Administra- 
tion (NTA) /Oslo  got  Network  Control  Center  (NCC)  permission 
to  break  the  communications  path  in  connection  with  replace- 
ment of  an  equalizer. 

The  SDAC  Communications  Circuit 

11/09  Loss  of  carrier  approx,  one  hour. 

12/05  NCC  suspected  line  trouble,  as  SDAC  claimed  dis- 

continuity in  data  transfer  between  the  two  data 
centers.  A modem  test  was  run  to  eliminate  possible 
malfunction  in  that  device. 

12/09  Trouble  with  the  system  most  of  the  day  from  0700  GMT. 

Otherwise,  fair  operation. 
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The  Terminal  Interface  Message  Processor  (TIP) 

TIP  preventive  maintenance  (PM)  carried  out  in  accordance 
with  Bolt,  Beranek  and  Newman  (BBN)  schedule.  With  a few 
exceptions  the  system  has  been  running  continuously  during 
most  of  the  period.  A Very  Distant  Host  (VDH)  test  was  imple- 
mented by  Mr.  Kelley  of  BBN  10  November  1977.  The  system  was 
taken  down  for  about  two  hours. 


We  have  also  experienced  teletype  trouble  a few  times.  5 January 
communication  with  NCC  was  impossible,  as  NCC  had  great  dif- 
ficulties in  reading  us.  Messages  addressed  to  NORSAR,  however, 
were  easily  interpreted.  The  interface  card  (in  the  CPU)  was 
the  possible  cause  for  the  trouble.  The  teletyne  was  nartly 
inoperative  in  February.  A card  was  replaced  (92/25)  . 

TIP  Connections 

No  change  in  the  connection  of  the  IMP  nart  of  the  TIP  since 
last  report.  On  the  other  hand,  a few  changes  have  been 
made  with  respect  to  TIP  port  (LIU)  connections.  Eleven  ports 
are  occupied,  of  which  the  Norv/egian  Defence  Research  Establish- 
ment (NDRE)  occupies  6 (Nos . 2,  4,  12,  40,  41  and  42),  NORSAR 
2 (Nos.  3 and  6),  Norwegian  Telegraph  Administration/Research 
Establishment  1 (No.  54),  University  of  Oslo  1 (No.  50),  Regne- 
anlegget  Blindern-K jeller  (RBK)  1 (No.  1). 


O . A . Hansen 
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III.  IMPROVEMENTS  AND  MODIFICATIONS 

111.1  Detection  Processor 

The  modifications  to  this  system  have  been  minor  within 
this  reporting  period.  The  system  seems  to  be  relatively 
bug-free,  and  the  two  errors  corrected  (see  below)  were 
not  critical  for  the  performance. 

As  was  discovered  at  SDAC,  the  subfield  identifications 
were  not  correct  in  the  result  record  from  the  Online 
Event  Processor,  transmitted  over  the  ARPANET.  This  was 
caused  by  a programming  error  in  the  PNRSAD  module, 
and  was  promptly  corrected  12  November. 

The  Alternate  Telemetry  Command  (ATC)  code  03  (Channel 
Gain  Measurement) , initiated  from  the  EOC,  gave  incorrect 
results.  This  was  traced  back  to  a missing  input  card 
for  the  Core  Image  Tape  Generation  program.  A re-run 
of  this  program,  with  the  missing  card  added,  produced 
a new  Core  Image  tape.  The  invalid  ATC  code  08  (and 
other  related  commands)  gave  proper  results  after 
system  restart  with  the  new  Core  Image  tape,  on  17 
November . 

111. 2 Event  Processor 

The  new  AUTOEP  processing  system,  which  reads  Online  Event 
Processor  (OEP)  results  off  the  Detection  Log  Tape  and  performs 
further  processing  on  these  data,  has  been  improved,  in  order 
to  meet  user  requirements.  The  program  may  now  be  stopped  grace- 
fully at  any  time,  giving  back  the  results  achieved  up  to  this 
point.  Also,  the  program  can  be  told  to  defer  processing  until 
OEP  results  later  than  a certain  time  have  been  read  in  from  the 
Detection  Log  tape.  Thus  the  analyst  may  stop  the  program  at 
will  and  restart  from  the  same  point  later.  A remaining  problem 
is  that  improper/invalid  data  infrequently  causes  a program 
check  in  the  filtering  routine.  Code  for  checking  traces  before 
filtering,  and  for  trapping  program  checks,  will,  however, 
shortly  be  implemented. 
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A version  of  this  program,  using  the  concurrent  plotting 
feature,  is  also  available.  When  this  version  is  running, 
event  plots  will  be  generated  at  the  same  tine  as  an  event 
is  processed.  In  cases  where  results  from  a sequence  of 
events  are  wanted  as  soon  as  possible,  this  feature  should  come 
in  handy. 

Since  manual  repunching  of  a large  proportion  of  the  produced 
bulletin  cards  from  this  program  seems  necessary,  work  to 
build  an  extension  to  this  system  has  started.  All  bulletin 
lines  produced  by  the  AUTOEP  program  are  now  automatically 
written  to  a Disk  Bulletin  file  (DBF)  at  the  end  of  processing. 

A new  program  will  access  the  DBF  and  display  bulletin  lines 
upon  the  screen  of  the  2260  Display  Station.  Lines  may  then 
be  modified  by  the  analyst.  An  important  new  feature  in 
this  program  will  be  the  automatic  recomputation  of  para- 
meters which  are  related  to  parameters  changed  by  the  analyst 
(i.e.,  change  in  arrival  time  gives  automatic  change  in  the 
origin  time).  A final  'publish'  will  produce  copies  of  the 
edited  bulletin  lines  on  cards  or  tape. 

D.  Rieber-flohn 

III. 3 Array  Instrumentation  and  Facilities 

A leftover  modification  from  1976,  namely,  modification  of 
BE-lightning  protection  cards  (Larsen  et  al,  1975)  was  com- 
pleted at  04C  in  October  1977. 

As  of  7 November  the  standard  low  pass  filters  with  upper 
3 dB  point  at  4.75  Hz  were  replaced  by  filters  with  upper 
3 dB  points  at  8.0  Hz  on  channels  01A06,  02B06,  02C06,  04C06 
and  06C06.  The  frequency  response  is  not  measured  exclusively 
for  these  channels,  for  practical  use  the  frequency  response 
of  NORSAR  Analog  SP  Station,  ref.  Fig.  III. 3.1  can  be  consulted. 
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From  15  March  a Teledyne  Geotech  S-500  seismometer  has  been 
in  operation  at  06C  channel  02,  uo  to  30  'larch  in  vertical 
position  and  thereafter  in  NS  horizontal  position.  The  con- 
nections are  given  in  Fig.  III. 3. 2.  An  internal  report  will 
be  issued  when  the  test  period  is  completed. 

A.  Kr . Nilsen 
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IV.  MAINTENANCE  ACTIVITY 

A brief  review  of  the  maintenance  accomnl ished  at  the  subarrays 
by  the  field  technicians  as  a result  of  the  remote  array 
monitorin'';  and  routine  inspections  is  given.  The  monitoring 
schedule  has  not  been  changed  in  the  oeriod. 

Maintenance  Visits 

Fig.  IV. 1 shows  the  number  of  visits  to  the  subarrays  in 
the  period.  The  subarrays  have  on  the  average  been  visited 
8.7  times  (without  01A,  4.5  times).  The  large  number  of  visits 
to  OlA  are  due  to  cable  breakages,  dryout  and  paintina  of  the 
LPV  and  LP  tanks  and  complete  checkout  of  the  LP  seismometers. 
Also  at  03 C the  LPV  and  LP  tanks  were  dried  out  and  painted. 

There  has  been  one  maintenance  visit  on  the  communications 
system  in  the  period. 


o 10  t i 


01A  01B  02B  02C  03C  04C  06C 

Fig.  IV. 1 Number  of  visits  to  the  NORSAR  subarrays  in  the  period 
1 October  1977  to  31  March  1978. 
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Preventive  Maintenance  Projects 

The  preventive  maintenance  work  in  the  array  is  described 
in  Table  IV. 1.  The  adjustments  are  corrections  of  character- 
istics within  the  tolerance  limits. 


Table  IV. 1 

Preventive  Maintenance  Work  in  the  Period 
1 October  1977  to  31  March  1978 
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Disclosed  Malfunctions  on  Instrumentation  and  Electronics 

Table  IV. 2 aives  the  number  of  accomplished  adjustments  and 
replacements  of  field  equipment  in  the  arrav  with  the  excep- 
tion of  those  mentioned  in  Table  IV. 1. 


Table  IV. 2 

Total  number  of  required  adjustments  and  replacements  in  the  NORSAR 
data  channels  and  SLEM  electronics 
(1  Oct  1977  - 31  Mar  1978) 


Malfunction  of  Rectifiers,  Power  Loss,  Cable  Breakages 


There  has  been  no  malfunction  of  the  rectifiers  in  the  period 
The  number  of  cable  breakages  was  two,  requiring  6 days'  work 
by  the  field  technicians. 

Array  Status 

Average  values  of  some  of  the  characteristics  are  given  in 
Tables  IV. 3 and  IV. 4,  compared  with  nominal  values. 

Table  IV.  3 

Average  values  of  the  channel  resolution  and 
channel  voltage  as  of  31  March  1978 


Subarray 

Channel  Resolution 

i SP 

L? 1 

PM/QU 

Voltage  P-P 

NM/QU 

Voltage  P-P 

01A 

42.14 

5.79 

2.48 

4.92 

01B 

41.27 

5.91 

2.54 

4.80 

02B 

41.80 

5.84 

2.61 

4.67 

02C 

[ 

40.94 

5.96 

2.61 

4.88 

03C 

40.80 

5.98 

2.53 

4.82 

04C 

42.10 

5.80 

2.56 

4.77 

06C 

45.52 

5.36 

2.60 

4.69 

Table  IV. 4 

Average  array  channel  characteristics  values 
as  of  31  March  1978 


Chan. 

Channel 

Resolution 

DC  Offset 

Nat.  Freq. 

Damping 

PM/QU 

Nominal 

Volts 

P-P 

Nomir  al 

Millivolts  Nominal 

Hz  Nominal 

Nominal 

SP 

42.09 

42.7 

5.80 

5.71 

-0.2  0 

1.03  1.00 

0.69  0.70 

LP 

NM/QU 

2 . 54 

2.47 

4.80 

4.94 

0 

20.0 

0.648  0.64 
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Conclusion 

Also  in  this  period  the  array  instrumentation  performance  has 
been  stable  and  satisfactory.  As  can  be  seen  under  the  section 
of  array  status  the  main  channel  characteristics  are  close 
to  nominal  and  with  little  change  from  previous  periods. 

Alf  Kr.  Nilsen 
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ABBREVIATIONS 


BB 

- 

Broad  band 

CAL 

- 

Calibration 

CMR 

- 

Common  mode  rejection 

CTV 

- 

Central  Terminal  Vault 

DC 

- 

Direct  current 

D CO 

- 

DC  offset 

EPU 

- 

External  power  unit 

FP 

- 

Free  period 

LP 

- 

Long  period 

LPV 

- 

Long  period  vault 

LTA 

- 

Line  Termination  Amplifier 

MP 

- 

Mass  position 

NM 

- 

Nanometer 

PM 

- 

Pi come ter 

P-P 

- 

Peak-to-peak 

QU 

- 

Quantum  units 

RA-5 

- 

SP  seismograph  amplifier 

RCD 

- 

Remote  centering  device 

SLEM 

Seismic  short  and  long  period  electronics 

module 

SP 

- 

Short  period 
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V.  DOC  UHL'.  NT  AT  I ON  DEVELOPED 

V . 1 Reports,  Papers 

Christof fersson.  A.,  1978:  Statistical  models  for  seismic 

magnitude,  NORSAR  Scientific  Report  1-77/78,  March  78. 
Gjtfystdal,  H.,  1977  : Final  Technical  Summary,  MOP.SAR  Scientific 
Rep.  No.  3-76/77,  Nov.  77. 

Haddon,  R.A.W.,  and  E.S.  Husebye,  1978:  Joint  interpretation 

of  P-wave  time  and  amplitude  anomalies  in  terms  of  litho- 
spheric heterogeneities,  Geophys.  J.R.  Astr.  Soc . , in  oress. 
Husebye,  E.S.,  P.C.  England  and  I.B.  Eanberq,  1978:  The  ideal- 
body  concept  in  interpretation  of  the  Oslo  Rift  gravity 
data  and  their  correlation  with  seismic  observations, 

In:  I.B.  Ramberg  and  E.  Neumann  (eds.):  The  Oslo  Paleo- 
rift,  NATO  ASI  Proceedings, 

Rieber-.Mohn , D.,  1978:  The  use  of  ARPANET  for  transmission  of 
real  time  seismic  data,  NORSAR  Internal  Rep.  No.  1-77/78, 
NTNF/NORSAR,  Kjeller,  Norway. 

Rieber-Mohn,  D. , 1978:  Documentation  of  the  NCP  task,  NORSAR 
Internal  Rep.  No.  2-77/78,  NTNF/NORSAR,  Kjeller,  Norway. 
Sacks,  I.S.,  A.  Snoke  and  E.S.  Husebye,  1R7R : Lithospheric 
thickness  beneath  the  Baltic  Shield,  Tec tonophysics , 
in  press. 

T;j0stheim,  D.,  1978:  Autoregressive  modelling  and  spectral 

analysis  of  array  data  in  the  plane,  Geophys.,  in  press. 
Tj0stheim,  D.,  and  O.A.  Sandvin,  1978:  Multivariate  auto- 
regressive feature  extraction  and  the  recognition  of 
multichannel  waveforms,  IEEE  Trans,  on  Comp.,  in  press. 

L.B.  Tronrud 

V.2  Program  Documentation 

Two  documents  have  been  completed  during  this  period,  describing 
how  plotting  can  be  performed  concurrently  with  the  plot- 
generating program. 

N/PD-91  describes  the  new  foreground  plotting  program  FDPLOT, 
which  cooperates  with  the  background  program  and  receives 
plot  data  from  it,  via  a shared  disk  file,  before  passinq 
these  to  the  plotter. 
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N/PD-92  describes  the  modified  PLOT  subroutine,  that  writes 
the  plot  data  to  the  shared  disk  file,  and  communicates  with 
the  foreground  (PLOT  program,  using  the  interpartition 
signalling  facility  in  the  operating  system. 

D.  Rieber-Mohn 
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V £ . SUMMARY  Of  SPECIAL  TECHNICAL  RE  PORTS/PAPERS  P7E^  REn 

In  this  section  a brief  summary  is  qiven  of  the  results  of 
ongoing  and  recently  completed  research  projects  at  NTtFVNORSAR. 
The  presentation  covers  research  conducted  under  NTNF/NORSAR' s 
contract  with  ARPA  as  well  as  research  nrojects  sponsored  by 
i Norwegian  authorities.  Of  oarticular  interest  to  the  seismic 

discrimination  problem  among  the  Norweeian-funded  undertakings 
is  the  participation  of  two  NTNF/NORSAP.  seismologists  in  the 
seismological  expert  qroup  established  by  the  United  Nations. 

In  addition,  some  of  the  research  conducted  at  NTN^/NORSAR 
in  connection  with  seismic  risk  studies  is  also  of  general 
seismological  interest,  and  is  therefore  included  in  the 
following . 


VI. 1 Work  of  the  Seismological  Expert  Croup  Established 
by  the  United  Nations 

On  22  July  1976  the  Conference  of  the  Committee  on  Disarmament 
(the  CCD)  of  the  United  Nations  established  an  Ad  Hoc  grouo 
of  Government-appointed  experts  to  consider  and  report  on 
international  cooperative  measures  to  detect  and  identify 
seismic  events,  so  as  to  facilitate  the  nonitoriro  of  a 
comprehensive  test  ban.  Representatives  of  a total  °7 
nations  participated  in  the  expert  group,  which  met  in 
Geneva,  Switzerland,  in  five  sessions.  Its  final  renort 
was  transmitted  to  the  CCD  on  9 March  19  73  and  contained, 
specific  recommendations  for  a global  system.  In  short,  the 
main  elements  of  the  recommended  system  were: 

(i)  A systematic  improvement  of  the  observations  reported 
from  a network  of  more  than  fifty  seismological 
observatories  around  the  globe. 

(ii)  An  international  exchange  of  these  data  over  the 
Global  Telecommunications  System  of  the  World 
Meteorological  Organization. 

(iii)  Processing  of  the  data  at  special  international 

data  centers  for  the  use  of  the  participant  states. 

The  report  also  considered  some  stems,  such  as  an  exnerimental 
exercise,  which  could  be  taken  initially  to  assist  the  establish- 
ment of  such  a cooperative  data  exchange  system. 
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The  Norwegian  government  appointed  Dr.  E.S.  Husebye  and 
Dr.  F.  Ringdal,  both  of  NTNF/NORSAR,  to  represent  Norwav  in 
the  expert  group.  Dr.  Ringdal  was  chosen  by  the  group  to  act 
as  its  scientific  secretary.  While  the  particination  in  the 
expert  group  was  funded  by  Norwegian  authorities,  oart  or  the 
research  work  done  at  NTNF/NORSAR  in  this  connection  has  also 
been  of  relevance  to  the  NTNF/NORSAR ' s ARPA  contract.  For 
example,  a comprehensive  detectability  study  o^  nearly  sgn 
globally  distributed  seismograph  stations  was  undertaken, 
and  the  results  have  now  been  published  (Rincdal  et  al, 

1977)  . 

For  supplementary  comments  on  the  work  of.  the  Ad  Hoc  rrrouo, 
we  refer  to  the  editorial  of  Nature,  6 April  lrPP  (see  Fig. 
VI. 1.1),  where  political  and  scientific  implications  of  the 
proposed  measures  are  discussed. 

E. S.  Husebye 

F.  Ringdal 
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6 April  1978 


Twenty  years  of  test  ban  talk 


Is  I9SX.  a conference  of  scientific  experts  in  Geneva 
made  tile  first  steps  lowaids  devising  an  international 
seismic  monitoring  system  vvlneli  would  verify  com- 
pliance with  any  ticatv  hanmirg  underground  nuclear 
weapons'  ’esls  In  |9Mt,  with  political  interest  in  a com- 
prehensive test  ban  in  the  doldmnis  but  with  ten  years 
of  scismological  research  v»n  a national  basis  completed, 
SIPR1,  the  Stockholm  International  Peace  Research 
Institute,  convened  further  informal  meetings  of 
scientists  m an  attempt  to  get  a comprehensive  test  ban 
((1 H)  talked  ahoul  again.  Now  in  1078.  with  serious 
political  discussion  proceeding  both  at  superpower  level 
and  amongst  a wide  range  of  nations  at  the  UN  Con- 
ference of  the  Committee  on  Disarmament  (CCD), 
scientists  have  again  reported  on  what  must  be  done  in 
an  intemation.il  contest  to  monitor  a test  ban  treaty 
1 lieu  report,  the  icsult  ol  deliberations  by  scientists 
from  27  countries  over  a period  of  a year  and  a half, 
b is  recently  been  released  bis  CCD  document  55S).  It 
reflects  substantial  credit  on  its  participant,  especially 
on  Dr  Ulf  Unc  -on  fiom  Sweden,  its  chairman,  and 
Dr  f Ringdal  from  Norway,  its  scientific  secretary,  for 
although  there  vva'  a clear  need  to  hammer  out  some 
form  of  consensus  in  the  document,  this  has  not 
prevented  its  message  from  ' being  clear  and  unam- 
biguous. 

The  science  of  test  ban  monitoring  was  mostly  done 
in  the  1960s.  rcchnitpies  to  increase  dctcctabilitv.  to 
discriminate  between  explosions  and  earthquakes,  to 
relate  seismic  magnitude  to  explosive  yield,  to  locate 
events  more  accurately  were  all  developed  rapidly  during 
that  period,  and  have  in  recent  years  undergone  rela- 
tively little  further  change.  W hat  hav  happened  in  the 
past  ten  years,  however,  ha-  been  a marked  improve- 
ment in  data  handling.  Studies  which  used  to  take 
months  of  data  accumulation  and  hand  measurement 
can  now  be  done  in  a morning  at  a computer  console 
Main  international  communication  links,  both  formal 
and  informal,  now  exist  ami  more  are  planned.  This,  of 
course,  is  true  in  many  other  blanches  of  science  and 
greatly  benefits  research,  hut  in  seismology  the  bonus 
is  lh.it  it  is  now  possible-  to  talk  of  an  international 
centre  or  ccillics.  vvuih  lapal  access  to  data  of  a high 
quality  from  seismomele's  all  round  the  world,  provid- 
ing a routine  (low  ol  information  highly  relevant  to 
the  verification  of  a ( lit  In  many  ways  the  recent 
repot t is  a hliicpiiiit  lor  such  ail  operation,  which  might 
he  pi  ended  by  an  experiment  taking  up  to  two  years. 

Ii  is  of  interest  to  compute  the  predictions  of  network 
capability  which  :u  ■ being  made  ill  197S  with  those  pul 


forward  in  I c>68  (which  came  essentially  from  a pro- 
digital  era).  The  detection  of  events  almost  invariably 
depends  on  the  successful  registration  above  noise  levels 
of  so-called  body  waves.  Hus  detection  capability  has 
improved  roughly  threefold;  explosions  of  yields  of  I or 
2 kilotons  in  hard  rock  in  most  parts  of  the  Northern 
Hemisphere  would  now  most  likely  be  picked  up.  The 
improvement  in  detection  of  surface  waves,  necessary 
to  the  identification  of  explosions  as  such,  is  even 
gi eater.  Identification  might  now  be  possible  for  shots 
as  low  as  5 to  10  kilotons  in  hard  rock. 

Not  all  the  progress,  however,  is  in  the  science  and 
technology.  For  the  past  twenty  years  the  Soviet  Union's 
willingness  to  co-operate  in  a scheme  of  test-ban  moni- 
toring has  been  in  doubt.  Many  times  she  h is  declared 
that  site  is  perfectly  prepared  to  sign  a treaty,  but  that 
she  regards  national  means'  as  adequate  for  verification 
Since  the  Soviet  national  seismic  network  is  of  very 
limited  value  in  monitoring  the  United  Slates,  this  state- 
ment is  open  to  the  interpretation  that  tiie  nature  of 
US  society  is  such  that  clandestine  small -yield  testing 
would  be  impossible.  But  the  corollary  is  that  the  nature 
of  Soviet  society,  and  even  the  geography,  leaves  the 
door  open  to  violation  and  that  much  wider  open  if  the 
Soviet  Union  will  provide  no  data  to  international 
agencies.  It  is  too  little  realised  that  at  present  even  the 
informal  channels  by  which  seismologists  exchange  data 
are  closed  on  the  days  that  the  Soviet  Union  conducts 
an  underground  lest. 

The  recent  discussions,  however,  offer  some  promise, 
flic  Soviet  Union,  a rather  hesitant  participant  to  begin 
with,  eventually  co-operated  fully,  and  even  allowed 
five  of  its  own  stations  to  be  used  in  various  calculations 
- in  contrast  to  the  French  and  Chinese  who  stayed 
away.  The  next  step  will  be  when  data  from  these  five 
vtations  are  supplied  on  a routine  basis.  This  is  unlikely 
to  happen  before  a treats  is  signed - the  Soviet  Union 
would  regard  provision  of  such  material,  containing 
possible  evidence  of  weapons  tests,  as  tantamount  to 
handing  out  state  secrets.  But  if  the  long-term  intention 
is  to  participate  fully,  this  must  be  regarded  as  an 
optimistic  sign 

A comprehensive  lest  han  needs  much  more  Ilian  a 
good  verification  network  to  hring  it  into  being.  But  this 
report  is  hound  to  provide  some  reassurance,  particu- 
larly in  tile-  United  States,  that  such  a network,  includ- 
ing Soviet  stations,  is  possible.  T lie  proposals  will  n-  t 
guai. mice  that  ls-sts  at  the  kiloton  level  can  be  positively 
identified  as  such.  But  they  do  show  some  cvidcncc- 
for  the  first  lime-  of  truly  international  goodwill  r] 
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VI . 2 Statistical  Models  for  Seismic  Magnitude 

The  concept  of  seismic  magnitude  - a measure  the  kinetic 
energy  of  the  elastic  waves  released  by  an  earthquake  ••  was 
first  suggested  by  C.F.  Richter  in  1936  (see  Richter,  195*?). 
Magnitude  measurements,  which  give  an  indication  of  the 
relative  size  of  earthquakes,  are  today  routinely  made  at  all 
seismological  observatories  and  represent  an  integral  part 
of  many  research  investigations.  In  the  context  of  seismic 
event  classification,  i.e.,  discrimination  between  earthquakes 
and  underground  nuclear  explosions,  the  magnitude  parameter 
is  of  paramount  importance.  The  reason  is  simoly  that  desoite 
extensive  research  efforts  the  so-called  m,  :M  discriminant 

D S 

is  still  considered  the  most  reliable  one  and  is  also  the 
most  widely  used.  On  the  other  hand,  in  certain  branches 
of  seismology  like  source  mechanism  studies  the  parameter 
seismic  moment  has  replaced  magnitude  (to  a large  extent) 
for  indicating  the  size  of  large  earthquakes  (e.q.,  see 
Kanamori  and  Anderson,  1975)  . In  other  parts  of  seismology 
the  magnitude  parameter  has  been  somewhat  discredited  because 
of  its  considerable  variability  due  to  different  physical 
factors,  some  of  which  are  difficult  or  impossible  to  quantify. 
It  must,  however,  be  said  that  some  of  the  research  aimed  at 
the  magnitude  problem  must  be  rated  as  rather  orinitive. 

In  view  of  the  importance  of  the  magnitude  parameter  in  a 
seismic  discrimination  context,  NORSAR  scientists  have  in 
recent  years  given  a considerable  attention  to  the  magnitude 
problem.  The  study  has  been  focused  on: 

(i)  whether  parts  of  the  observed  magnitude  scattering 
was  associated  with  inhomogeneities  in  the  earth 
(forward  scattering  of  small-scale  inhomogeneities) , 

(ii)  the  magnitude  estimation  itself,  and  finally 

(iii)  developing  discriminants  having  a better  performance 
than  that  of  the  m^ifl 
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This  section,  together  v/ith  section  VI.  3,  describes  sore 
recent  efforts  regarding  the  first  two  factors  mentioned 
above,  namely,  the  scattering  of  the  P-wave  amplitude 
observations  and  in  particular  the  proper  estimation  of 
m^-magni tudes  given  the  observations  from  a network  of  seis- 
mograph stations  and  arrays.  Extensive  studies  shot'  that 
the  P-wave  amplitude  variations  across  the  hORSAR  array  are 
rather  large  and  may  be  clearly  associated  with  structural 
heterogeneities  at  the  bottom  of  the  lithosphere  (see  Sec. 

VI. 3).  This  amplitude  scattering  has  a relatively  short 
wavelength,  i.e.,  varyinc  rapidly  with  small  changes  in  dis- 
tance and  azimuth.  Also,  the  amplitude  distribution  across 
the  NORSAR  array  may  be  approximated  by  a loanormal  statistical 
distribution.  This  behavior  has  also  been  observed  for  world- 
wide amplitude  data,  and  implies  that  the  station  naanitude 
correction  term  and  thus  the  scattering  term  in  maanitude 
estimation  models  can  be  considered  a Gaussian  variable. 

A novel  approach  to  the  magnitude  estimation  problem  was 
the  work  of  Ringdal  (1976),  who  introduced  a maximum  likelihood 
technique  for  estimating  magnitude  from  a network  of  stations, 
thereby  taking  into  account  information  on  stations  being 
operational,  but  not  detectino  weaker  events.  Ignoring  the 
latter  kind  of  information  would  in  most  cases  result  in  a 
positive  magnitude  bias  for  small  events. 

The  mentioned  maximum  likelihood  approach  has  recently  been 
further  elaborated  by  Christof f ersson  (1978).  His  approach 
differs  from  that  of  Ringdal  (1976)  in  that  it  takes  into 
account  the  probability  that  the  event  is  actually  detected 
by  the  network,  whereas  Ringdal  (1976)  considered,  in  statistical 
terms,  a sample  space  which  also  included  cases  where  an 
event  was  not  seen  by  any  of  the  stations  in  the  network. 

The  practical  difference  in  the  estimates  provided  by  the 
two  methods  is  generally  small,  and  the  relative  merit  of 
the  two  approaches  will  not  be  discussed  here. 
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In  summary,  the  statistical  models  presented  in  the  oaper 
by  Christof fersson  (1973)  in  connection  with  seismic  magnitude 
deals  with  two  main  situations.  The  first  concerns  the  estima- 
tion of  magnitude  for  an  event  using  a fixed  network  of 
stations  and  taking  into  account  the  detection  and  bias 
properties  of  the  individual  stations.  The  second  treats  the 
problem  of  estimating  seismicity  and  detection  and  bias  proper- 
ties of  individual  stations.  The  models  are  aoolied  to  analyse 
the  magnitude  bias  effects  for  an  earthquake  aftershock 
sequence  from  Japan,  as  recorded  by  a hypothetical  network 
of  15  stations.  It  is  found  that  network  magnitudes  computed 
by  the  conventional  averaging  technique  are  considerably 
biased,  and  that  a maximum  likelihood  approach  using  instan- 
taneous noise  level  estimates  for  non-detecting  stations 
gives  the  most  consistent  magnitude  estimates.  Finally, 
the  models  are  applied  to  evaluate  the  detection  character- 
istics and  associated  seismicity  as  recorded  by  three  vela 
arrays  (UBO,  TFO,  WHO) . 

While  the  two  statistical  situations  discussed  by  Christof fersson 
(1978)  each  provide  powerful  techniques  for  eliminating  the 
bias  caused  by  non-detections  of  individual  stations  of  a 
network,  they  are  in  general  suited  for  two  different  estima- 
tion problems.  The  first  (or  conditional)  approach  is  useful 
mainly  for  estimating  the  magnitude  of  individual  events,  and 
can  be  applied  equally  well  to  earthquakes  and  explosions. 

The  second  (or  unconditional)  approach,  gives  a convenient 
framework  for  joint  estimation  of  structural  parameters  such 
as  seismicity  (a  and  b in  the  recurrence  formula  log  N=a-b*?!) 
and  station  bias.  It  can  also  be  used  to  estimate  station 
detection  characteristics.  This  second  method  is  a unified 
approach  which  provides  a generalization  of  earlier  works 
of  Kelly  and  Lacoss  (1969)  and  Ringdal  (1975) . 
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Perspectives . The  various  maximum  likelihood  approaches 
discussed  above  for  ensuring  consistent  maanitude  estimates 
have  to  our  knowledge  only  been  applied  to  m^-observations . 

The  main  reason  for  this  is  that  this  kind  of  maanitude  data 
are  the  only  ones  which  are  easily  and  abundantly  available. 

There  is  no  reason  why  these  novel  estimation  techniques 
should  not  be  applicable/extended  to  surface  wave  maanitude 
(M  ) estimation,  and  also  other  sinnificant  problems  in  a 
discrimination  context  like  the  m^:M  relationship  in  particular 
for  weak  events.  Research  on  these  types  of  nroblems  is  now 
in  progress,  and  our  efforts  here  are  concentrated  on  one 
hand  on  developing  algorithms  where,  for  examnle,  oossible 
correlation  between  the  P and  Rayleigh  wave  detectability 
for  a given  station  is  taken  into  account,  and  on  the  other 
hand  to  constructing  comprehensive  data  bases  from 

both  array  and  SRO-recordings . The  use  of  advanced  statistical 
techniques  in  analyzing  the  relationship:  is  expected 

to  give  more  definite  answers  to  a number  of  questions  as 
to  the  nature  of  the  relationship  e.g.,  the  ranae  over  which 
it  may  be  considered  linear,  the  associated  slope  (both  for 
explosion  and  earthquake  populations)  and  most  importantly, 
its  behavior  at  low  magnitudes  where  the  nroblems  due  to  non- 
detections are  considered  to  be  most  significant. 

A.  Christof f ersson , Univ.  of  Upnsala 

F.  Ringdal 

E.S.  Husebye 
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VI. 3 Short  Period  P-wave  Amplitude  Variability 

In  an  earlier  NORSAR  technical  summary  report  (Gjpystdal,  1977) 
we  described  in  some  detail  a rather  extensive  experiment  on 
a deterministic  approach  to  the  analysis  of  short  period  P-wave 
amplitude  anomalies  as  observed  across  the  NORSA.R  array.  The 
essence  of  this  work  was  that  from  observed  travel  time  residuals 
various  models  for  lithospheric  velocity  contrasts  in  the 
form  of  thin  lenses  were  constructed.  Then  usinc  a finite  dif- 
ference calculation  scheme  the  expected  P-wave  amplitude 
distribution  was  obtained  which  in  turn  was  compared  to  that 
actually  observed  from  real  events.  The  best  fit  between  ob- 
served and  calculated  amplitude  anomalies  was  obtained  for 
lenses  in  the  depth  range  150-200  km,  although  the  theoretical 
amplitude  values  could  only  account  for  about  some  40%  of  the 
variance  in  the  observational  data.  However,  the  fit  between 
observed  and  theoretical  results  is  clearly  better  than  the 
quoted  number  indicates,  that  is,  the  respective  anomaly 
patterns  are  quite  similar  and  the  indicated  deteriorated  fit 
stems  partly  from  problems  in  'exact'  positioning  of  the  two 
patterns.  A reasonable  physical  explanation  here  appears  to 
be  that  the  observed  wavelength  of  observed  P-time  anomalies 
is  significantly  larger  than  that  of  corresponding  P-amplitude 
anomalies  which  implies  that  the  ' time ' -derived  lens  models 

are  not  sufficiently  detailed  for  very  precise  amplitude  < 

calculation.  A consequence  of  this  hypothesis  is  that  the 

reverse  process  should  be  more  rewarding,  that  is,  we  would 

intuitively  expect  a better  fit  if  we  tried  to  predict  time 

anomalies.  This  has  actually  be  achieved  using  an  energy 

flux  formulation  for  the  lens  focusing/defocusing  effects 

which  ultimately  led  to  a Poisson  differential  equation.  The 

main  results  here  were  a correlation  of  about  0.90  between 

I 

observed  and  predicted  time  anomalies.  i 

I 


On  the  basis  of  this  study  (Haddon  and  Husebye,  1978)  it  is 
concluded  that  time  and  amplitude  anomalies  originate  from  the 
same  lithospheric  structures  which,  as  a good  first  approximation, 
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may  be  represented  in  terms  of  a 2-D  heterogeneous  layer  at 
depths  around  150-200  km  or  the  bottom  of  the  lithosphere. 

Vie  note  in  passing  that  the  lithosphere,  an  intecrral  concent 
of  modern  plate  tectonics,  is  not  well  defined  seismolooica] 1 v . 
However,  more  recent  observations,  also  at  NORSAR,  of  S-to-P 
converted  waves  indicated  a well-defined  discontinuity  at 
a depth  of  around  230  km  (Sacks  et  al,  1973J  see  also  Sec. 

VI.  11)  which  may  be  taken  to  indicate  a lithospheric  thickness 
of  the  same  order.  Furthermore,  velocity  oerturbations  re- 
quired for  accounting  for  the  anomalous  P-wave  amplitude 
observations  amounts  to  a few  per  cent  and  thus  are  directly 
compatible  with  similar  results  obtained  by  Aki  et  al  (1976, 
1977)  and  by  scattering  analysis  of  precursor  and  coda 
waves  (Husebye  et  al,  1976,’  King  et  al , 1975,'  Haddon  et  al, 

1977)  . 

Perspectives . The  above  study  has  been  completed  (Haddon  and 
Husebye,  1978)  and  some  wider  applications  have  been  considered, 
that  is,  can  the  methodology  used  here  be  adapted  to  other 
arrays,  networks  or  conventional  seismograph  stations.  Indeed, 
we  have  undertaken  some  preliminary  analysis  of  LASA  data 
but  as  the  spatial  earthquake  sampling  of  this  array  is  less 
symmetrical  than  that  of  NORSAR  we  have  temporarily  halted 
this  v/ork.  As  pointed  out  above,  a major  problem  in  reproducing 
amplitude  anomalies  from  time  anomalies  is  the  difference  in 
wavelengths  of  these  two  types  of  anomalies.  On  the  other 
hand,  the  time  residual  projection  scheme  used  in  constructing 
the  lens  models  has  proved  very  valuable  in  analysis  of 
absolute  travel  time  anomalies  (say  those  listed  in  ISC- 
catalogues)  . For  example,  using  absolute  NOP.SAR  time  anomalies 
we  have  reproduced  those  areas  of  overlap  of  the  thin  lens 
models  derived  from  relative  travel  times.  This  result  was 
not  too  unexpected  in  view  of  the  exceptionally  hiah  quality 
of  the  NORSAR  data,  but  as  demonstrated  by  Husebye  and 
Ringdal  (1978),  the  above  projection  scheme  may  also  be  an 
alternative  to  conventional  analysis  of  time  residuals  from 
seismograph  networks  of  continental  dimensions.  Furthermore, 
we  are  also  considering  a joint  inversion  scheme  of  NORSAR  time 


and  amplitude  data  and  in  th„s  oarticular  case  based  cn  ray 
tracing  principles.  Finally,  we  do  consider  this  tyne  or 
problem,  that  is,  a better  unders'.andino  oF  intrinsic  ampli- 
tude variations  in  particular  for  near-field  observations  in 
the  distance  range  0°-30°  to  be  of  fundamental  imnortance  i 
a seismic  discrimination  context.  A key  word  to  success  here 
is  of  course  flexibility  in  modellirq  seismic  heterogeneities 
in  the  mantle. 

E.S.  Hu-ebyc 

R.A.VJ.  Haddor.  (L’niv.  of  Sydney) 
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VI. 4 Inversion  of  Travel  Time  Data 

The  NORSAR  interest  in  geophysical  inversion  problems  dates 
back  to  the  summer  of  1974  when  Professor  K.  Aki,  MIT, 
visited  Kjeller.  The  research  then  initiated  resulted  in 
development  of  a novel  and  in  particular  flexible  technique 
for  inversion  of  travel  time  residuals  so  as  to  produce  a 
3-D  image  of  the  seismic  velocity  anomalies  beneath  the  array 
or  network  in  question  (e.g.,  see  Aki,  Christof f ersson 
and  Husebye  (ACH) , 1976,  1977,  and  Husebve  et  al,  1976). 

This  particular  inversion  technique  has  become  rather  popular 
for  detailed  studies  of  lithospheric  heterogeneities  in 
various  parts  of  the  U.S.  and  also  Japan  and  even  adapted 
to  inversion  of  time  residuals  from  the  global  seismographic 
network  (Dziewonski  et  al,  1977). 

In  the  amplitude  inversion  experiment  described  in  Sec. 

VI . 3 we  mentioned  that  the  construction  of  the  thin  lenses 
used  in  calculating  theoretical  amplitude  values  was  based 
on  a relatively  simple  projection  scheme  which  resulted  in 
a 2-D  seismic  velocity  anomaly  model.  Furthermore,  Haddon 
and  Husebye  (1978)  used  essentially  the  same  data  as  Aki  et  al 
(1977)  (see  VI. 3)  so  apparently  the  two  mentioned  studies 
gavj  conflicting  results  and/or  the  bulk  of  lithospheric 
inhomogeneities  are  confined  to  a relatively  thin  layer 
in  the  lower  part  of  the  lithosphere.  We  do  consider  that  the 
differences  between  the  Aki  et  al  (1977)  and  Haddon  and  Husebye 
(1978)  studies  can  be  partly  reconciled  by  using  more  homo- 
geneous model  specifications  and  partly  reflect  a fundamental 
problem  in  seismology,  namely,  that  of  discriminating  between 
a relatively  thick,  weakly  inhomogeneous  layer  and  on  the 
other  hand  a relatively  thin,  strongly  inhomogeneous  layer. 
Furthermore,  in  the  3-D  inversion  scheme  the  basic  model 
parameters  like  number  and  thicknesses  of  layers,  average 
layer  velocities  and  block  sizes  are  not  subject  to  estima- 
tion but  are  specified.  For  example,  the  ACH-inversion 
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scheme  can  easily  reproduce  the  Haddon  and  Husebye  lower 
lithospheric  lenses  using  a 2-layer  model  with  the  second 
one  located  at  depths  around  150-200  km,  and  still  have 
roughly  the  same  variance  reduction  as  obtained  within  the 
original  5- layer  model  used  in  the  ACH-publication . Indeed, 
the  difference  between  the  ACH-results  and  the  Haddon- Husebye 
lens  models  is  relatively  minor  as  the  3 bottom  layers  in 
the  ACH-model  have  very  similar  velocity  anomaly  patterns  which 
in  turn  are  very  similar  to  that  of  the  lens  models.  Besides 
the  possibility  that  the  velocity  anomalies  in  the  lithosphere 
in  the  NORSAR  siting  area  may  have  a significant  vertical 
extent,  the  specification  of  the  ACH  basic  model  parameters 
may  have  some  important  bearing  on  the  final  results  and 
consequently  on  their  subsequent  interpretation.  Part  oc  the 
problems  here  are  intuitively  obvious  as  the  standard  errors 
of  the  estimated  velocity  anomalies  are  relatively  larger 
thus  implying  that  the  physical  resolution  may  be  less  than 
generally  assumed.  This  problem  was  indeed  discussed  when 
the  ACH-paper  was  written,  but  at  that  time  hampered  by 
limited  accessability  to  sufficiently  larqe  computers  for 
running  the  inversion  program. 

In  view  of  the  apparent  controversy  between  the  ACH-results 
and  those  of  Haddon  and  Husebye,  and  also  the  popularity 
of  the  ACH- inversion  technique,  we  have  initiated  research 
in  this  problem  and  have  so  far  designed  a scheme  by 
which  v/e  can  simulate  a large  class  of  basic  model  specifica- 
tion parameters  by  only  solving  a limited  number  of  the 
total  number  of  linear  equations  involved.  So  far  the  cor- 
responding programming  efforts  have  been  minor,  but  the 
computer  programs  are  expected  to  be  completed  in  the  near 
future. 

E.S.  Husebye 

A.  Christof fersson  (Univ.  of 
Uppsala) 
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VI. 5 Microearthquake  Surveillance  of  Svalbard 

Up  until  now  only  very  limited  date,  have  beer  available 
for  the  study  of  the  seismic  activity  in  and  around  Svalbard, 
and  most  of  the  studies  so  far  have  been  based  on  teJeseisnic 
data  (Hodgson  et  al,  1965;  Sykes,  1965  J Husebye  et.  al,  1975," 
Bungum,  1977).  However,  even  i^  only  about  1-2  intraylate 
earthquakes  from  Svalbard  are  reported  every  year  from  the 
teieseismic  recordings,  the  large  (M  =5.9)  event  in  the 
Storfjorden  area  on  13  January  1976  clearly  showed  that  a 
certain  seismic  hazard  is  present.  The  fact  that  this  event 
had  a faulting  mechanism  atypical  for  intraolate  earthquakes 
(Bungum,  1977)  also  emphasized  the  need  for  a close1-  investi- 
gation, which  could  be  done  only  by  installina  seismic  stations 
on  the  archipelago  itself. 

A program  for  such  microearthquake  surveillance  o^  the  Svalbard 
Archipelago  v/as  initiated  by  installation  of  seismic  stations 
in  Barentsburg  (BBG)  , Longyearbyen  (LYR)  and  Pyramiden  (P^D)  , 
in  cooperation  between  the  Russian  mining  trust  1 Arktikuool ' , 
Store  IJorske  Spitsbergen  Kullkompani,  the  Norwegian  Polar 
Institute  and  NTNF/NORSAR  (Bungum  et  al,  1978).  The  installation 
of  the  seismometers  (Sprengnether  MEQ-300)  was  done  in 
December  1977,  and  the  operation  during  the  first  few  weeks 
v/as  somewhat  unstable,"  this  was  due  to  various  technical 
problems  most  of  them  instabilities  related  to  installation 
in  sub-zero  temperatures.  For  this  reason,  reliable  time 
corrections  were  not  available  for  the  data  analyzed  in  this 
report.  When  analyzing  the  first  2k  months  of  data,  we 
consequently  had  to  use  only  the  S--P  times,  from  which  locations 
were  calculated  using  a maximum  likelihood  procedure  which 
also  uses  all  available  information  about  the  various  errors 
involved.  (For  details,  see  Section  VI. 7).  The  adooted 
relation  between  the  S-P  times  and  epicentral  distance  v/as 
based  upon  the  results  of  Mitchell  et  al , 1979. 
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Seismograms  have  been  available  and  read  from  the  three 
microearthquakc  stations  for  various  time  intervals  between 
8 December  1977  and  24  February  1978  and  from  KBS  for 
December  and  January.  A total  of  687  earthquakes  have  been 
detected  at  one  or  more  of  the  stations,  515  of  which  (or 
75%)  are  local  events.  The  best  station  is  where  an 

average  of  7.5  local  events  per  day  (corrected  for  down  times) 
have  been  detected.  The  number  for  PhD  is  4.3  and  for  BBO  ^ . 8 , 
whereas  the  poorest  station  in  this  respect  is  the  SSN 
station  KBS,  where  only  2.0  local  events  per  day  were  detected 
(see  also  Section  VI. 6) . Magnitudes  have  been  computed  for 
231  of  the  events,  and  they  all  fall  in  the  ran^e  0-3,  with 
a peak  at  around  magnitude  1.0. 

The  total  number  of  located  events  is  234.  About  one  half  o^ 
the  events  are  located  using  2 stations  and  the  other  hal^ 
with  3 stations.  The  locations  of  the  latter  ones  are  shown 
as  an  epicenter  map  in  Fig.  VI.  5.1,  where  a <-"-eat  cluster 
, of  events  at  the  west  side  of  Storfjorden  appears  as  the 

dominating  feature.  The  precisions  of  the  locations  are  so 
far  not  good  enough  for  a closer  delineation  of  this  highly 
active  earthquake  zone,  since  the  size  of  the  cluster  is 
not  much  larger  than  the  computed  uncertainty  elliDse  for 
each  event. 

The  eastern  coast  of  Nest  Spitsbergen  in  Storfjorden  thus 
appears  to  be  an  area  of  high  intraplate  seismic  activity, 
since  an  average  of  7.5  events  per  day  with  a probable 
origin  in  Storfjorden  have  been  recorded  by  the  instrument 
in  Longyearbyen.  The  work  of  Mitchell  et  al  (197'’)  showed 
. most  of  the  epicenters  to  be  confined  within  a narrow  E-W 

trending  zone  at  77.7°N  about  30  km  long  and  suggested  that 
' , the  fault  plane  of  the  18  January  1976  earthquake  (M  =5.9) 

in  this  area  (Bungum,  1977)  was  along  this  seismic  zone. 

The  state  of  stress  in  the  crust  and  delineation  of  active 
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zones  of  weakness  where  faulting  may  take  place  are  valuable 
information  for  present  and  future  industrial  development 
on  Spitsbergen.  Spitsbergen  is  cut  by  major  fault  zones 
dating  back  to  Paleozoic  with  no  associated  teleseismically 
recorded  activity  (Husebye  et  al,  1975) . A future  semi- 
permanent network  of  microearthquake  stations  on  Spitsbergen 
would  be  capable  of  detecting  zones  that  are  tectonically 
active  but  have  a low  level  of  seismic  activity. 

A further  step  towards  this  end  will  be  taken  by  installa- 
tion of  a new  station  (with  digital  recording)  in  Svea 
(S4  in  Fig.  VI. 5.1)  in  May  1978,  while  we  hope  to  be  able 
to  install  two  more  stations  at  a later  stage  (S5  and  S6 
in  Fig . VI . 5 . 1)  . 

H.  Bungum 

H.  Gjtfystdal 

B.  Kr.  Hokland 

Y.  Kristof  fersen  , 'Jorv/ectian 
Polar  Institute 
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Fig.  VI. 5.1 


Epicenter  map  of  earthquakes  located  on  the  basis  of 
the  three  stations  BBG,  PRD  and  LYR.  The  future  station 
in  Sveagruva  is  also  indicated  (S4)  , as  well  as  two 
possible  sites  on  the  west  side  of  Storfjorden  (S5  and 
S6) , that  if  available  would  greatly  improve  the  array 
configuration . 
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VI . 6 Teleseismic  Detectability  of  the  Svalbard  Microearthquake 
Network 

As  discussed  in  Section  VI . 5 , the  primary  puroose  of  the  seis- 
mograph network  installed  on  Svalbard  north  of  Norway  has  been  to 
survey  the  local  seismic  activity.  Nonetheless,  we  have  also 
investigated  the  teleseismic  performance  of  each  station,  and 
preliminary  results  are  encouraginr.  Table  vi . 6 . 1 lists  the 
geographical  coordinates  of  the  four  seismoaraoh  stations 
presently  in  operation  on  Svalbard,  and  Table  VI. 6.°  shows 
their  relative  detection  performance  over  a 2^-month 
period.  Not  unexpectedly,  each  of  the  three  microearthauake 
stations  detect  considerably  more  events  than  the  IT7SSN 
station  at  Kings  Bay  (KBS) . More  interesting  is  that  this 
appears  to  hold  true  not  only  for  local  earthquakes,  but 
also  with  regard  to  teleseismic  events.  (Averacres  of  2.2, 

2.1  and  1.5  teleseismic  events/day  versus  1.4  at  KBS). 

Because  of  the  short  time  interval  available,  and  the  lack 
of  complete  overlap  of  the  periods  under  consideration, 
these  numbers  should  be  considered  only  to  give  tentative 
indications.  We  plan  to  conduct  a more  comprehensive  detect- 
ability study  at  a later  stage,  usinq  more  complete  data. 
Meanwhile,  we  note  that  the  high  sensitivity  of  the  micro- 
earthquake seismographs  to  high  frequency  signals  make  these 
instruments  particularly  suitable  to  detect  underground 
explosions.  As  an  example.  Fig.  VI. 6.1  shows  the  P-wave  recorded 
at  station  PRD  from  a presumed  explosion  in  Eastern  Kazakh, 

26  March  1973,  with  m^ (NORSAR) =5 . 2 . Another  Eastern  Kazakh 
presumed  explosion  from  19  March  (m^  (NOP.SAP)  =5 . 1)  was  also 
detected,  although  the  signal-to-noise  ratio  was  considerably 
lower  than  in  the  former  case.  The  epicentral  distance  to 
Eastern  Kazakh  is  about  40  degrees  ’,  in  comparison,  the  dis- 
tance from  Svalbard  to  Novaya  Zemlya  is  only  10  dearees. 

We  plan  to  expand  the  study  of  teleseismic  and  near-field 
detectability  of  the  Svalbard  network  as  more  data  become 
available. 

F.  Pingdal 

H.  Bunaum 

B.  Kr.  Hokland 
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Tab! e VI.C.l 

Names  and  coordinates  of  the  seismic  stations  used  in  this  study.  In  the 
last  column  is  given  the  percentage  of  the  time  between  8 December  1977 
and  24  February  1978,  during  which  each  of  the  stations  has  been  in  ooera- 
tion  (or  data  available  in  case  of  KBS) . 


Site 

Code 

Lat 

Long 

Data 

Availabilty  (%) 

Operated  by 

Barentsburg 

BBG 

78.073 

14.240 

89.5 

Norwegian  Polar 

Pyramiden 

PRD 

78.659 

16.303 

58.7 

Institute  & 
NTNF/NORSAR 

Longyearbyen 

LYR 

78.189 

15.578 

63.8 

Kings  Bay 

KBS 

78.918 

11.924 

69.6 

Univ.  of  Bergen 

i 


Table  VI. 6. 2 

Detectability  statistics  for  the  four  stations  used  in  this  study.  Data 
for  the  stations  BBG,  PRD  and  LYR  cover  the  time  period  between  8 December 
1977  and  24  February  1978,  whereas  data  from  KBS  have  been  available 


only  for  December  1977  and  January 
corrected  for  station  down-time. 

1978. 

The  daily 

averages  have 

been 

BBG 

PRD 

LYR 

KBS 

Total 

Detected  events,  total 

449 

303 

482 

185 

687 

- Average,  per  day 

6.3 

6.5 

9.6 

3.4 

Detected  events,  local 

336 

199 

376 

111 

514 

- Average  per  day 

4.8 

4.3 

7.5 

2.0 

Detected  events,  teleseismic 

113 

104 

106 

74 

173 

- Average  per  day 

1.5 

2.2 

2.1 

1.4 

I 


I 


Best 

Available 

Copy 
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Fig.  VI. 6.1  P-wave  recording  at  PRD,  Svalbard,  of  a presumed 

explosion  from  Eastern  Kazakh,  26  March  1978.  Signal 
onset  is  shown  by  an  arrow. 
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VI.  7 A Maximum  Likelihood  Procedure  for  Local  Event  Location 
Based  on  Observed  S-P  Time  Differences  at  Two  or  More 
Stations 

A method  has  been  developed  in  order  to  compute  local  event 
locations  merely  based  on  the  relative  arrival  times  of  P 
and  S waves  observed  at  a number  of  individual  stations. 

The  method,  which  has  been  applied  on  data  from  the  Svalbard 
microearthquake  network  (see  section  VI. 5,  and  Bunqum  et  al, 
1978)  takes  advantage  of  the  fact  that  there  normally  is 
a constant  ratio  between  P and  S velocities  in  the  crust, 
making  (for  short  distances)  the  epicentral  distance  A 
approximately  linearly  dependent  upon  the  S-P  time: 

A =*  k • t(S-P)  (1) 

Knowing  the  epicentral  distance  from  two  stations,  we  may 
usually  compute  two  epicenters  symmetrically  located  about 
the  line  connecting  the  stations.  Having  a distance  observation 
from  one  or  more  additional  stations  located  non-symmetrically 
relative  to  the  former  ones,  we  will  generally  be  able  to  choose 
the  proper  solution,  however,  in  this  case  the  final  location 
should  be  based  on  a sort  of  'averaging  process'  since  the 
'distance  circles'  will  normally  not  intersect  each  other 
in  one  single  point,  due  to  the  distance  errors  involved. 

The  present  location  procedure  is  based  on  the  maximum 
liklihood  principle  from  statistical  theory.  Assuming  that 
the  error  in  the  'observed'  epicentral  distance  A^  for  a 
given  station  i is  normally  distributed  with  zero  mean  and 
standard  deviation  cm  , we  may  locally  (close  to  the  epicenter) 
approximate  the  'distance  circle'  by  a straight  line  and 
express  the  associated  probability  density  function  as  a 
'Gaussian  ridge'  distributed  about  this  line  (see  ^ig. 

VI. 7.1)  : 
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Here,  x and  y are  rectangular  coordinates  centered  in  a point 
in  the  vicinity  of  the  true  epicenter,  and  a^,  and  c^  are 
parameters  defining  the  'distance  line'  in  this  coordinate 
system. 


Having  chosen  the  origin  of  this  system,  f.ex.,  in  the  inter- 
section point  between  two  arbitrarily  chosen  distance  circles, 
the  parameters  a^,  b^  and  c^  may  be  easily  computed  from 
the  station  coordinates  and  the  'observed'  value  of  the  distance 

V 

When  an  expression  like  (2)  has  been  found  for  N stations, 
we  may  compute  the  joint  probability  density  of  the  epicenter 
by  forming  the  product 


N 

p(x,y)  = H p. (x , y ) (3) 

i=l  1 

and  locate  the  epicenter  in  the  point  corresponding  to  the 
maximum  value  of  p(x,y)  which  can  be  shown  to  represent  a 
binormal  distribution  for  N>2.  In  addition  to  the  location 
of  the  maximum  point  (point  of  maximum  likelihood) , we  can 
analytically  find  the  axes  and  orientation  of  the  confidence 
ellipses  of  the  resulting  distribution. 

H.  Gjpystdal 
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Fig.  VI. 7.1 


Error  distribution  of  the  ’distance  line’  from  each 
station,  calculated  from  the  P-S  times. 
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VI. 8 Precisely  Located  Earthquakes  in  the  Vicinity  of  NO^SAR 

Even  though  local  earthquakes  in  the  siting  area  o^  NONPAR 
are  quite  rare  (Husebye  et  al,  197P) , there  occurred  between 
1970  and  1978  four  such  events  that  could  be  subjected  to 
a detailed  hypocentral  investigation  using  the  NOR8AR  record- 
ings. 

The  locations  were  performed  using  the  program  HYP071  (Lee 
and  Lahr,  1975) , and  the  essential  nart  of  the  experiment 
was  concentrated  on  the  estimation  of  a useful  crustal 
model.  For  this  purpose,  three  local  explosions  for  which 
the  locations  were  precisely  known  were  used,  and  we  thereby 
came  up  with  the  following  model : 

Depth  (km)  Velocity  (km/s) 


0-16 

6.25 

16-32 

6.65 

>32 

8.15 

Using  this  model  and  the  above-mentioned  location  orogram,  j 

the  hypocenters  given  in  Table  VI. 8.1  were  estimated.  The  I 

estimated  standard  epicentral  error  for  all  four  events  was 
around  or  less  than  one  kilometer. 

t 

I 

| 

Table  VI. 8.1  ' 

Estimated  hypocentral  coordinates  for  four  earthquakes  in  the  NORSAR  area.  i 


No. 

Date 

Time  of  Day 

Lat  (N) 

Long  (E) 

Depth 

1 

71.07.19 

00.59.11.5 

60°43 . 1 1 

10°43 .6 ' 

31 

2 

73.10.01 

16.44.14.3 

59°55. 1 ' 

11°25.4' 

6 

3 

73.11.23 

06.49.36.9 

60°33 . 2 1 

11°28 . 2 ' 

23 

4 

77.12.11 

21.46.12.1 

60°56. 7 ' 

10°53 . 1 ' 

22 

6 

23 

22 
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in  order  to  i nvestigatc  the  resolution  of  the  deoth  estima- 
tion, an  experiment  was  performed  in  which  the  hyoocenter 
was  constrained  at  given  depths  between  0 and  4n  kn  and  the 
associated  RMS  error  in  seconds  was  computed.  The  results 
are  shown  in  Fig.  VI. 8.1,  where  it  is  seen  that  the  resolu- 
tion is  fairly  good  for  all  of  the  events. 

The  hypocentral  solutions  presented  here  are  the  most  accurate 
ones  so  far  published  for  earthquakes  in  pennoscandia , 
and  this  applies  in  particular  to  the  focal  denths.  It  is 
therefore  interesting  to  note  that  our  results  are  consistent 
with  those  of  Husebye  et  al  (1978),  who  found  from  macroseismic 
data  that  most  of  the  depths  should  be  in  the  ranee  15-30  km. 

H . Bunqum 
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VI. 9 Hacroseismic  Data  Collection  Using  Newspaper  Ads 


On  the  basis  of  the  promising  results  obtained  in  developing 
a method  and  a procedure  for  automatic  processing  of  macro- 
seismic  questionnaires  (Husebye  et  al,  1976),  it  was  found 
desirable  to  change  the  format  of  the  questionnaires  so  that 
they  would  be  more  directly  suited  for  analysis.  Also,  with 
the  present  great  interest  shown  by  the  news  media  in  connec- 
tion with  local  earthquakes,  it  was  also  tempting  to  take 
advantage  of  this  professionally.  Consequently,  a project 
was  initiated  (in  cooperation  with  University  of  Bergen) 
in  which  macroseismic  questionnaires  now  are  regularly 
published  as  newspaper  ads,  in  a format  shown  in  Tig.  VI. 9.1. 
The  questions  are  answered  by  crossing  for  'yes',  'no'  or 
'don't  know'.  So  far,  the  new  format  and  the  new  data 
collection  procedure  have  been  tested  on  the  four  earthquakes 
listed  in  Table  VI. 9.1,  where  it  is  seen  that  230  replies 
were  received  for  one  of  the  events.  While  the  data  still  are 
under  analysis  for  all  of  the  events.  Fig.  VI. 9. 2 shows  the 
geographical  distribution  of  the  answers,  which  of  course  in 
this  case  is  convolved  with  the  distribution  of  the  news- 
papers . 


H.  Bungum 
E.S.  Husebye 
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Table  VI. 9.1 

List  of  earthquakes  for  which  questionnaires  have  been  published  as 

newspaper  ads . 


No.  Date  Place  (in  Norway) 

1 11  Dec  1977  Brumunddal 

2 9 Jan  1978  Stord 

3 20  Mar  1978  M0re 


No.  of  Replies 

108 

69 

230 


29  Apr  1978 


Grong 


Jordskjelvet  ved  M0re 
den  20.  mars  kl.  04.58 

Som  et  ledd  l et  forskningsprosjekt  i snmarbeid  mellom  det  selsmologiske  observatoriet  NOR- 
SAR  p4  Kjeller  or  Jordskjelvstasionen  1 Bergen  ber  dlsse  lnstitusjonene  her  om  pubUkums 
hjelp  U1  4 kartlegge  virknlngene  av  jordskjelvet  den  20.  mars  kl.  04.58.  - Sporsm41ene  retter  seg 
til  beboerne  av  dot  enkelte  hus/Iellighet,  or  dersom  sporsm41et  lkke  er  aktuelt,  krysses  det  av  1 
kolonnen  for  -vet  ikke-.  En  er  ogsA  interessert  1 svar  fra  personer  som  bor  1 rtmelig  naerhet  av 
skjelvet  og  som  ikke  bar  merket  noe. 

Navn  og  adresse  (bare  hvis  onskrsl  

Noyaktig  stedsanglvel.se  (vlktig!):  

Jordbunnstype.  Sand  □ beire  □ bosmasser  (lkke  speslflsert)  □ Kompakt  grunn  (sand, 
stein)  □ Fjell  □ 

Bygnlngstype:  Tre  □ Mur/lett  hetong  □ Betong  ['] 

Ja  Nel  Vet  lkke 

□ □ □ Ble  skjelvet  merket’  (Hvis  • Nel-  oiler  -.Vet  ikke»,  kan  ovrlge  sporsm41  Ignoreres). 

□ □ □ Ble  skjelvet  merket  av  de  fleste  1 husct/leiligheten? 

□ □ □ Ble  skjelvet  merket  av  alle  1 buset/lelllgheten? 

□ □ □ Ble  skjelvet  merket  utendors? 

□ □ □ Ble  skjelvet  merket  av  de  fleste  utendors? 

□ □ □ Ble  rystelsen  merket.  som  en  svak  skjelving? 

□ □ □ Ble  vibrasjonene  merket  som  en  passerende  lastebil/tog? 

□ □ □ Ble  rystelsen  merket  som  ved  en  sprengnlng? 

□ □ □ Hvis  om  natten,  ble  sovende  personer  vekket? 

□ □ □ Var  det  vanskelig  4 holde  balansen  mens  skjelvet  p4glkk? 

ODD  Kllrret  vlndusruter  eller  ovner? 

□ □ □ Svinget  tamper  eller  andre  hengende  gjenstander? 

□ □ □ Knaket  det  1 vegger  eller  gulv? 

□ □ □ Ristet  mobler? 

□ □ □ Beveget  bilder  p4  veegen  seg? 

□ □ □ Stoppet  pendelur? 

□ □ □ Ble  smft  og  lettc  gjenstander  forflyttet  (eller  veltet)? 

□ □ □ Var  det  noen  dorer  eller  vinduer  som  4pnet  seg  eller  slo  lgjen? 

□ □ □ Fait  boker  ned  fra  hyllnne? 

□ □ □ Ble  storre  mobler  forflyttet? 

□ □ □ Ble  glass  eller  porselen  knust’ 
n T]  [1  Ble  mobler  beskadlget? 

li  □ f ] Ble  smft  sknder  pft  hus  eller  murvegger  observert? 

□ U (.]  Ble  det  observert  storre  .sknder  pA  hus  eller  murvegger? 

U □ □ Ble  det  observer!  skader  pA  voler  (sprekker  I velbancn)? 

□ □ □ Ble  det  observert  sknder  pA  vannlednlngcr? 

Vennllgst  kryss  av,  kllpp  ut  og  returner  skjemaet  til: 

JORDSKJELVSTASJONISN,  Allcgt.  41,  5014  Bergen  U. 
eller:  NTNF/NORSAR,  Posthoks  51,  2007  Kjeller. 


Fig.  VI. 9.1  An  example  of  an  earthquake  questionnaire  as  a newspaper 
ad,  20  March  1978. 
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!><•■  i ;;rn  i <:  i ty  of  East  Africa 

As  part  of  a seismic  risk  analysis  for  a olanned  dan  oroiect 
in  Stiegler's  Gorge,  Tanzania,  NTNF/NORSAR  has  conducted 
an  extensive  study  of  the  seismicity  and  tectonics  of  Eastern 
Africa.  Seismologically , this  is  a very  interesting  region 
as  it  encompasses  the  East  African  Rift  System,  and  it  has 
consequently  been  the  subject  of  numerous  studies  in  the 
past  (see,  e.g.,  Maasha  and  Molnar,  197?).  Fia.  VI. 10.1 
shows  the  distribution  of  seismicity  in  the  area  covered 
by  the  present  investigation,  based  on  a catalogue  of  4n6R 
known  earthquakes  compiled  by  NTNF/NORSAR.  The  hich  seismic 
activity  along  the  rift  zones  may  be  clearly  identified) 
however,  it  is  of  interest  to  note  that  the  seismic  activitv 
shows  wide  distribution  also  outside  this  main  system. 

In  this  respect,  the  seismicity  of  East  Africa  shows  clear 
similarities  to  what  has  been  observed  in  other  intranlate 
areas  such  as  Northern  Europe,  Russia  and  Northeast  America. 

Special  attention  was  given  in  the  NTNF/NOPSAR  studv  to 
comparing  the  earthquake  magnitudes  reported  by  various 
agencies  to  those  calculated  from  NORSAR  recordinos.  Fig. 

VI. 10. 2 shows  a plot  of  NORSAR  versus  PDE  (Preliminary 
Determination  of  Epicenters,  U.S.  Oeolonical  Survey)  reported 
magnitudes.  Although  the  data  base  is  limited  due  to  the 
short  time  period  covered  (6  years) , it  is  evident  that, 
relatively  speaking,  PDE  shows  a systematic  positive  magni- 
tude bias,  which  is  most  pronounced  for  NORSAR  magnitudes  of 
about  4.0.  We  attribute  this  to  the  network  maanitude  bias 
problem  discussed,  e.g.,  by  Ringdal  (1976).  Similar  results 
were  found  when  comparing  NORSAR  m^  to  those  of  other  anencies 
such  as  the  International  Seismological  Centre  (ISC) . This 
is  consistent  with  observations  from  other  reaions,  and  it 

* 0 

appears  that  earthquakes  with  a PDE  or  ISC  renorted  m^  of 
5.0  or  above  in  many  cases  will  have  a considerable  oositive 
. ' bias,  sometimes  as  much  as  a full  magnitude  unit  relative  to 

a hypothetical  'true'  magnitude. 

F.  F.inndal 
H.  Bunoum 
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Fig.  VI. 10.1  Locations  of  all  4060  events  in  the  NTNF/NORSAR 
earthquake  catalogue  for  Eastern  Africa. 
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VI. 11  Lithospheric  Thicknesses  in  the  General  NORSAR 
Siting  Area 

An  integral  concept  in  modern  plate  tectonic  hypothesis 
is  the  structural  units  lithosphere  and  asthenosphere, 
and  as  such  also  widely  used  in  seismological  contexts. 

The  basic  differences  betwen  these  two  uppermost  layers  of 
the  earth  are  slightly  lower  velocity  but  much  lower  attenua- 
tion and  viscosity  in  the  asthenosphere  than  in  the  litho- 
sphere. The  a-factor  is  of  particular  importance  in  near- 
field (range  0-30°)  seismological  studies,  as  event  detect- 
ability would  be  inversely  proporitional  to  the  Q-factor. 
Furthermore,  there  are  considerable  regional  lithosphere 
differences;  for  example,  heat  flow  and  surface  wave  studies 
indicate  that  the  lithosphere  in  shield  areas  like  Fenno- 
scandia  is  relatively  thick  as  compared  to  oceanic  areas 
(Pollack  and  Chapman,  1977,’  Lee  and  Solomon,  1975) 

1974).  In  the  latter  case,  the  lithosphere-asthenosphere 
transition  is  generally  marked  with  a pronounced  velocity 
reversal.  The  combination  of  thick  lithosphere  and  high 
Q-values  clearly  implies  good  seismic  event  detectability 
in  the  near-field  distance  range  as  demonstrated  by 
Khalturin  et  al  (1977)  and  some  preliminary  results  on 
related  problems  have  been  published  by  Husebye  et  al 
(1977) . Anyway,  the  topic  of  this  section  is  to  describe 
an  experiment  aimed  at  estimating  lithospheric  thickness 
by  observations  and  analysis  of  so-called  S-to-P  converted 
waves  associated  with  a discontinuity  tentatively  interpreted 
as  the  lithosphere/asthenosphere  boundary  (see  Fig.  VI. 11.1). 
This  research  started  in  1976  when  Dr.  I.S.  Sacks  of 
Carnegie  Institute,  Washington,  D.C.,  visited  NOPSAF,  and 
has  now  been  completed  (see  Sacks  et  al,  1978). 


Arrivals  interpreted  as  Sjd  with  a conversion  beneath  the 
Baltic  Shield  (Umea  area,  NE  Sweden)  were  found  at  NORSAR 
for  5 events  in  the  epicentral  ranqe  70-82°  with  back 
azimuths  of  37-57°.  In  the  following  we  discuss  the  analysis 
of  one  of  these  events  in  some  detail  with  emphasis  on 
particle  motion  processing  of  the  records  (see  Husebye  et 
al,  1975).  Discontinuities  in  the  particle  motion  can  be  used 
to  determine  the  initial  onset  of  phases  whose  lonq  period 
character  makes  this  difficult  to  establish  precisely  from 
the  time  domain  record.  Fig.  VI. 11. 2 shows  the  vertical 
radial  and  transverse  components  and  the  particle  motion 
for  the  minute  preceding  direct  S.  Before  Sjd  (Fig.  VI. 11. 2 (1)) 
the  particle  motion  is  not  in  the  earthquake-station  great 
circle  path.  After  the  Sjd  onset  the  particle  motion  is  along 
azimuth  (Fig.  VI. 11. 2 (3  upper)  and  in  the  vertical  plane 
the  particle  motion  (Fig.  VI. 11. 2 (3  lower))  is  a tiaht  el- 
lipse. The  onset  of  Sjd  can  be  determined  by  tracing  the  motion 
backwards  in  time  to  the  point  where  the  particle  motion  breaks 
from  the  smooth  elliptical  motion  (Fig.  VI. 11. 2 (2  lower)).  The 
motion  associated  with  Sjd  persists  until  the  shear  arrival 
(Fig.  VI. 11. 2 (4)),  when  the  dominantly  radial  motion  changes 
to  transverse. 

Using  arrival  times  determined  from  the  particle  motion,  the 
differential  travel  time  for  this  event  is  found  to  be 
29+1  sec  which  implies  a cnversion  depth  of  250  + 15  km. 

Very  similar  results  were  obtained  for  the  other  4 events. 

I. S.  Sacks  (Carnegie  Inst.) 

E.S.  Husebye 

J. A.  Snoke  (Virginia  Poly- 
technical  Institute  & 

State  University) 
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